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The CDMS low ionization threshold experiment (CDMSlite) uses cryogenic germanium detectors 
operated at a relatively high bias voltage to amplify the phonon signal in the search for weakly 
interacting massive particles (WIMPs). Results are presented from the second CDMSlite run with 
an exposure of 70 kg days, which reached an energy threshold for electron recoils as low as 56 eV. 

A fiducialization cut reduces backgrounds below those previously reported by CDMSlite. New 
parameter space for the WIMP-nucleon spin-independent cross section is excluded for WIMP masses 
between 1.6 and 5.5 GeV/c®. 

PACS numbers: 95.35.+d, 14.80.Ly, 29.40.Wk, 95.55.Vj 


Gosmological and astrophysical measurements indicate 
that over one quarter of the energy density of the universe 
consists of nonbaryonic and nonluminous matter [1, 2]. 
Weakly interacting massive particles (WIMPs) produced 
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in the Big Bang are a compelling class of candidate par¬ 
ticles for this dark matter [3]. Recent accelerator [4] and 
direct-detection results [5] constrain the simplest super- 
symmetric models, sparking interest in alternative theo¬ 
ries, including theories with low-mass dark matter [6-11]. 
It has been suggested that, because the energy densities 
are similar for baryonic and dark matter, the relic den¬ 
sity may be generated by an asymmetry related to the 
baryon asymmetry [12-16]. In this case, the number den¬ 
sities of the two are also related, suggesting searches for 
particles with masses of a few GeV/c^. Hints of WIMP 
signals near detector thresholds [17-21], and an excess 
of gamma-ray emission from the Galactic Genter [22], 
have also contributed to interest in the low-mass WIMP 
region [5, 23-32]. 

WIMPs may create keV-scale recoils in laboratory de¬ 
tectors by elastically scattering off target nuclei [33]. 
WIMPs significantly below the target nuclei mass are an 
experimental challenge because they deposit small recoil 
energies, making it harder to distinguish WIMP signals 
from background and electronic noise. At these energies, 
techniques used to discriminate between nuclear recoils 
(NRs) and electron recoils (ERs) often diminish in ef¬ 
fectiveness. The GDMS low ionization threshold exper¬ 
iment (CDMSlite) [23] exchanges NR discrimination for 
a lower threshold. This paper describes data taken with 
one SuperGDMS iZIP detector [34] operated in GDMSlite 
mode at the Soudan Underground Laboratory [35]. 

GDMSlite employs the Neganov-Luke effect [36-38], 
which amplifies phonon signals when electric charges drift 
through a material [39-41]. Work done by the electric 
field on the electron-hole pairs is fully converted into 
phonons along the drift path of the charges. The to¬ 
tal measured phonon signal E^, including the energy E^. 
from the initial particle interaction, is 

Et = Er + Neh^Vb , (1) 

where e is the elementary charge, 14 is the bias voltage, 
and iVeh is the number of electron-hole pairs created in 
the interaction. The average energy required to create 
an electron-hole pair for an electron recoil in germanium 
is e.y = 3 eV/pair, which gives iVeh = Erje^. When I 4 is 
large compared to e^, the Neganov-Luke phonons domi¬ 
nate the measured signal and allow for lower thresholds 
to be reached. 

Nuclear recoils produce electron-hole pairs less effi¬ 
ciently, so an ionization yield Y(Er) that depends on 
energy and interaction type is dehned through iVeh = 
Y{Er)Er/e^, where Y =1 for electron recoils. The total 
energy can be expressed as 

Et = Er(^l + YiEr)^y (2) 

The energy scale thus also depends on the interaction 
type. The detector is calibrated using ERs and the re¬ 
sulting energy scale is labeled as keVee. The spectrum 
is then converted to nuclear recoil equivalent energy (la¬ 
beled as keVnr) by comparing Eq. 2 for ERs and NRs 


and solving for E^,-: 


1 + eVb/e 


1 + Y{En,)eVble^ 


( 3 ) 


where Y[E^r) is the ionization yield for NRs. 

During a first short data run of GDMSlite, this 
mode of operation was proven viable and provided lead¬ 
ing sensitivity to WIMPs with masses between 3 and 
6 GeV/c^ [23]. The data presented here are from the 
second CDMSlite run using the same detector, performed 
from February to November 2014 and taken in three hour 
long data series (to maintain the neutralization of the 
detector [35]). A voltage bias of —70 V was applied to 
one side of the detector with the other side at 0 V. The 
electronics setup followed that of the first CDMSlite run 
with phonon sensors instrumented only on the grounded 
side of the detector [23]. The single CDMSlite detec¬ 
tor was still part of the full SuperGDMS detector array 
where all detectors are read-out with every trigger. Sev¬ 
eral hardware and operational improvements were im¬ 
plemented for this run [42]. Fluctuations in bias volt¬ 
age were reduced by cleaning the high-voltage biasing- 
electronics board, sealing it against humidity, and plac¬ 
ing it in a dry nitrogen atmosphere. The large detector 
leakage current observed at the start of each data series in 
the first run was decreased by “prebiasing” the detector 
at —80 V for 10 min prior to each series. Cryocooler- 
induced microphonic noise that limited the threshold of 
the first CDMSlite run was better rejected by the in¬ 
stallation of vibration sensors near the connection of the 
cryocooler to the SuperGDMS cryostat. 

The cryocooler degraded during the run, causing the 
induced noise to dominate the trigger rate. The exper¬ 
iment was warmed to room temperature in July to al¬ 
low for the routine replacement of the cryocooler cold 
head. The run resumed in September and, due to the 
maintenance reducing the microphonic noise rate, a lower 
energy threshold was achieved. The run is thus natu¬ 
rally split into two periods: February-July (period 1) 
and September-November (period 2). 

Time intervals with exceptionally high trigger rates 
were removed from the WIMP-search exposure. Events 
with elevated prepulse noise were conservatively removed 
because the integrity of the detector cannot be guaran¬ 
teed during the time since the preceding event. Glitches, 
defined as pulses with uncharacteristically sharp rise and 
fall times, were observed and removed in three categories: 
electronic glitch events that cause simultaneous triggers 
in multiple detectors, glitches in the outer charge channel 
of the GDMSlite detector, which do not cause triggers 
themselves but can be coincident with a phonon trig¬ 
ger, and glitches that mimic standard events in all but 
pulse shape. Events coincident with the NuMI neutrino 
beam [43], and events whose NuMI timing information 
(for determining coincidence) was unavailable, were re¬ 
moved. The combination of the above cuts, except the 
third class of glitches (see below), reduced the live time 
from 132.23 to 115.59 d: 97.81 d (110.28 raw) in period 
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1 and 17.78 d (21.95 raw) in period 2. Events with en¬ 
ergy deposited in multiple detectors or coincident with 
the muon veto were also removed, with a combined effi¬ 
ciency of 97.21 ± 0.01 % for the detection of dark matter. 

Neutron calibrations were performed three times over 
the course of the run by exposing the detectors to a ^®^Cf 
source. Neutron capture on ’^°Ge creates ^^Ge, which de¬ 
cays via electron capture with a half-life of 11.43 d [44]. 
X rays and Auger electrons are emitted with a total en¬ 
ergy corresponding to the ^^Ga electron binding energy 
of the shell from which the electron is captured. The K- 
, L-, and M-shell binding energies are 10.37, 1.30, and 
0.16 keV, respectively [45], with the latter two in the 
energy region of interest for this analysis. 

The ER energy scale was calibrated using the K- 
capture line from ^^Ge. Drifts in detector bias and cryo¬ 
stat base temperature caused this line to vary by 5-10% 
over time. After correcting for those effects, two small 
(^2.5 % total) residual shifts of unknown origin were ob¬ 
served and corrected. 

The trigger efficiency was calculated by measuring, as a 
function of the reconstructed energy in the GDMSlite de¬ 
tector, the fraction of events triggering another detector 
that also trigger the GDMSlite detector. The data used 
for this were a subset of the ^®^Gf calibration events, 
with very strict cuts removing events that could cause 
triggers due to noise or cross talk in the GDMSlite detec¬ 
tor in coincidence with the other iZlP detectors. In the 
two periods, 50 % trigger efficiency was reached at 75'^'^ 
and 56^4 eVee, respectively, except for the 4.22 days im¬ 
mediately prior to the midrun warm-up. During these 
few days, the hardware trigger threshold was increased 
slightly to reduce noise-induced triggers (see below). 

Pulse shape was used to distinguish signal events from 
noise events, to make fiducial cuts, and to correct en¬ 
ergy calibrations. All events were fit to three different 
templates corresponding to the standard signal event, 
electronic glitches with quick rise and fall times, and 
cryocooler-induced low-frequency noise (LF noise). The 
differences between goodness-of-fit quantities for the tem¬ 
plates were used to identify glitches and LF noise. As 
the cryocooler degraded during period 1, the rate of 
LF noise greatly increased. A metric based upon the 
noise profile across the cryocooler’s 830 ms cycle was 
used to identify periods of calendar time with high or 
low cryocooler-induced LF noise. The pulse-shape dis¬ 
crimination cut was then set to be tighter during pe¬ 
riods of more cryocooler noise, with the goal of min¬ 
imizing leakage at the expense of efficiency. The effi¬ 
ciency of all three pulse-shape criteria was determined 
by a Monte Garlo pulse simulation that combines exper¬ 
imental noise (taken throughout the entire run) with the 
standard pulse template, scaled to the desired energies. 
Figure 1 gives the result for all three pulse-shape cuts, 
with the dominant loss in efficiency caused by the LF 
noise cut. Because the energy-estimating algorithm as¬ 
signs energies of ~80 eVee to LF noise, the goodness-of- 
fit separation between LF noise and good pulses becomes 



FIG. 1. (color online) Binned total signal efficiency after se¬ 
quential application of selection criteria: single-scatter and 
muon-veto (orange dashed-dot curve), pulse-shape (blue dot¬ 
ted curve), hardware-trigger (green dashed curve), and radial 
fiducial-volume (black solid curve) criteria. The gray band 
around the final curve shows the combined uncertainty on 
the overall efficiency at Icr confidence. 

less distinct near this point. This requires a harder cut at 
this energy, leading to a sharp drop in efficiency. A small 
systematic uncertainty corresponding to the variation of 
pulse shape with energy and position was estimated by 
measuring the efficiency for a range of pulse-shape tem¬ 
plates. 

A valid pulse shape has two components: a fast one 
whose amplitude depends on the position of the scatter¬ 
ing event and a slow one that carries the primary energy 
information. For this analysis, a new algorithm was in¬ 
troduced that fits pulses from each phonon channel with 
a linear combination of fast and slow template pulses, 
allowing the position and energy information to be sep¬ 
arated. The position information was used, along with 
the segmentation of the phonon sensor into one outer and 
three inner segments, to construct a radial parameter by 
comparing the start time and the amplitude of the fast 
component of the outer and inner sensors. This radial 
parameter is shown as a function of energy in Fig. 2, 
where the densely populated band at higher parameter 
values corresponds to events in the outer part of the de¬ 
tector. A nonuniform field in this region draws charge 
carriers to the sidewall of the detector, preventing them 
from traversing the full potential. This produces a re¬ 
duced Neganov-Luke amplification and distorts the en¬ 
ergy spectrum. 

The radial cut removes more than 90 % of these events, 
along with a small contribution of low-energy surface 
events originating on the detector housing, down to low 
energies while maintaining a reasonable fraction of the 
exposure for inner events. The few remaining reduced- 
energy events contribute to the background at lower en¬ 
ergy. In period 2, a cluster of background events ap¬ 
peared below ~250 eVee and was located in the outer part 
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FIG. 2. Radial parameter as a function of energy for the first 
(top) and second (bottom) periods. The dashed lines indicate 
the radial cut. The densely populated band at larger values 
corresponds to events near the edge of the detector. The ver¬ 
tical clusters are the ^^Ge capture peaks. The clear separation 
between outer and inner events decreases at high energy due 
to signal saturation in the outer phonon channel. The slight 
downward shift in the distribution after the maintenance pe¬ 
riod (caused by a small change in the operating point of the 
phonon sensors) together with the appearance of an unex¬ 
plained localized background below ~250 eVee near the edge 
of the detector motivated a tighter radial cut in period 2. 

of the detector near one of the connectors. This, together 
with differences in the operating conditions between the 
two periods, motivated a tighter cut in period 2. The 
fiducial-volume cut significantly reduced the background 
rate compared to the first CDMSlite run. 

The acceptance for the radial fiducial cut was de¬ 
termined using the ^^Ge electron-capture events, which 
sample a uniform spatial distribution in the detector. 
These events can be divided into three categories: those 
degraded in energy, those with full energy that fail the 
fiducial cut, and those with full energy that pass the fidu¬ 
cial cut. The fraction of events in the first category is 
given by the electric-held geometry and is energy inde¬ 
pendent. To measure this effect, the radius-energy plane 
was divided into sections and a likelihood-based Monte 
Carlo simulation was applied to each section indepen¬ 
dently to determine the contribution of two components: 
a time-independent background and a contribution from 
the ^^Ge activation lines exponentially decaying in time. 
The known ratio of L- to X-capture rates was used to 
separate the L- and if-capture contributions. The frac¬ 
tion of events with a full Neganov-Luke phonon signal 
was determined to be ~86 %. 

Next, the fraction of events with full phonon signal 
removed by the radial cut was computed at the capture- 
peak locations as the number of events passing the cut 
criterion divided by all peak events after background sub¬ 
traction. The background in the inner part of the detec¬ 
tor is negligible compared to the peak rate; in the outer 
part the background was calculated from the observed 
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FIG. 3. Spectrum of events passing all selection criteria, cor¬ 
rected for all efficiencies except the trigger efficiency. Dashed 
lines indicate the prominent features of the ^^Ge electron- 
capture decay with peaks at 10.37 keV (IF-shell), 1.30 keV 
(L-shell) and 0.16 keV (M-shell). Inset: Englargement of 
the lowest energies that determine the low-mass WIMP sen¬ 
sitivity, including the L- and M-shell activation peaks. Solid 
vertical lines show the 50 % trigger-efficiency points for the 
two periods. 

event rates above and below the peak. To measure the 
efficiency at lower energies, a pulse-simulation method 
was implemented. All events from the L-capture peak 
(chosen to avoid observed signal saturation in the outer 
phonon channel above ^2 keVee) were used to generate 
nearly noise-free pulses using the extracted composition 
of the fast and slow templates. These noise-free pulses 
were then scaled to the desired energy before adding mea¬ 
sured noise. This sample of artificial raw events was an¬ 
alyzed in the same manner as the real raw data. The 
efficiency was measured using the fraction of artihcial 
events passing the radial cut, taking into account the 
background contribution in the original event sample. 
The combined fiducial-volume efficiency was calculated 
to be ~50 % with a mild energy dependence as shown in 
Fig. 1. 

The hnal spectrum after application of all selection 
criteria and correcting for all efficiencies (except the trig¬ 
ger efficiency) is shown in Fig. 3. The main features 
are the ^^Ge electron-capture peaks at 10.37, 1.30, and 
0.16 keV. Hints of other peaks can be seen on top of a 
smooth background from Gompton scattering of higher- 
energy gamma rays. The observed ratio of the rates for 
the M- and T-capture peaks is 0.16 ± 0.03, compared to 
an expected 0.17 [46]. Other numerical characterizations 
of the primary components of the spectrum are listed in 
Table I. 

In GDMSlite mode, the ionization yield cannot be mea¬ 
sured on an event-by-event basis, necessitating a model. 
The most common model in the field is that of Lind- 
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Energy 

[keVee] 

Resolution 

k/M,%] 

Range 

[keVee] 

Average rate 
[keVee kg day] 

0.16 

1.30 

10.37 

11.4±2.8 

2.36±0.15 

0.974±0.009 

0.056-0.14 

0.2-1.2 

1.4-10 

11-20 

16.33l?-^| 

1.09±0.18 

1.00±0.06 

0.30±0.03 


TABLE I. Left: resolution of the ^^Ge capture peaks. Right: 
average rate between the peaks, after application of all selec¬ 
tion criteria, corrected for efficiency. The difference in rates 
above and below the TL-capture peak can be attributed to un¬ 
resolved peaks due to cosmogenic backgrounds and the higher 
rate below the M-capture peak can be attributed to more 
background leakage at lower energies in period 2. 

hard [47]: 



Y{Er,r) 


k- g{e) 
l + k-g{e) ’ 


( 4 ) 


where g{e) = 3e°-^®-l-0.7e° ®-l-e, e = 11.5ii^nr(keV)Z“^/^, 
and Z is the atomic number of the material. For germa¬ 
nium, k = 0.157, but the model is somewhat uncertain 
for low recoil energies. This uncertainty was accounted 
for by varying k uniformly between k = 0.1 and 0.2, 
which encompasses a majority of the experimentally ob¬ 
served data [48]. 

The data below 2 keVee and above the respective 
50% trigger-efficiency values were used to set the 90% 
confidence upper limit on the spin-independent WIMP- 
nucleon cross section. This was done using the optimum 
interval method with no background subtraction [49], 
the Helm form factor, and the following standard dark 
matter halo assumptions: a local dark matter density 
of 0.3 GeVcm“^, a most probable Galactic WIMP ve¬ 
locity of 220 kms“^, a mean orbital velocity of Earth 
with respect to the Galactic Genter of 232 kms“^, and a 
Galactic escape velocity of 544 kms“^ [50]. The results 
are shown in Fig. 4. 

Statistical and systematic uncertainties were propa¬ 
gated into the limit by calculating the final efficiency 
as a function of energy numerous times, each time pick¬ 
ing at random from the distributions of each input pa¬ 
rameter. Statistical uncertainties exist in the trigger 
efficiency, pulse-shape and radial-cut simulations, and 
energy-independent cuts, while systematic uncertainties 
exist in the pulse-shape and radial-cut simulations, and 
the Lindhard model fc-value. Limits were computed us¬ 
ing 1000 sample curves with the median and 95 % inter¬ 
val given in Fig. 4. The uncertainty is dominated by the 
Lindhard model, particularly below masses of 3 GeV/c^, 
and the radial-cut efficiency. 

This result excludes new parameter space for WIMP 
masses between 1.6 and 5.5 GeV/c^. The improvement 
in sensitivity over the first GDMSlite run is due to the 
increase in exposure, the reduction in threshold, and the 
decrease in background resulting from the radial fiducial- 
volume cut. A kink in the limit is seen at ^6 GeV/c^. 
Simulations indicate that this feature is a consequence of 


FIG. 4. (color online) Median (90% C.L.) and 95% inter¬ 
val of the WIMP limit from this analysis (black thick solid 
surrounded by salmon-shaded band) compared to other se¬ 
lected results. Other 90 % upper limits shown are from the 
first GDMSlite run (red thin solid curve) [23], SuperCDMS 
(red thin dashed curve) [24], EDELWEISS-II (red thin dot¬ 
ted curve) [25], LUX (dark-yellow thick dashed-dot curve) [5], 
GRESST (magenta thick dashed curve) [27], and DAMIC 
(purple thick dotted curve) [28] . Glosed regions are CDMS II 
Si 90% G.L. (blue dashed shaded region) [17], and CoGeNT 
90% C.L. (dark-green shaded region) [19]. 

the M-shell line at 160 eVee- Finally, the effect of hav¬ 
ing a tighter radial threshold in period 2 was considered. 
Placing the same looser threshold in both periods would 
result in a ^9% weakened sensitivity below masses of 
6 GeV/c^, which is well within the presented uncertainty 
band. 

In conclusion, the second GDMSlite run was success¬ 
ful in operating an iZIP detector at a bias potential of 
—70 V for a 70.10 kg days analysis exposure, with ioniza¬ 
tion thresholds of 75 and 56 eVee attained for the first 
and second period, respectively. The development of a 
fiducial-volume cut reduced the overall background rate 
significantly. The results presented here can be signif¬ 
icantly improved in future GDMSlite runs by lowering 
the threshold and background rate. The former can be 
achieved with better phonon resolution and higher bias 
potentials, and the latter with material selection and 
quality control. All of these improvements are planned 
for the SuperCDMS SNOLAB experiment. 
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